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Abstract 

We study the magnetic properties of a mixed Ising ferrimagnetic 
system, in which the two interacting sublattices have spins cr, (±1/2) 
and spins S, (±3/2, ±1/2) in the presence of a random crystal field, 
with the mean field approach. The obtained results show the exis- 
tence of some interesting phenomena, such as the appearance of a 
new ferrimagnetic phase namely the partly ferrimagnetic phase {m^ = 
^ , ms = +1) and consequently the existence of three topologically dif- 
ferent types of phase diagrams. The effect of increasing the exchange 
interaction parameter J, at very low temperature is investigated. The 
transitions shown in these phase diagrams are in good agreement with 
those obtained in the ground state case. 

PACS: 05.50.+q; 75.10Hk;75.50.Gg; 
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1 Introduction 

The magnetic properties of two sublattices mixed spins —1/2 and spin 
S > 1/2 Ising system, with a crystal field interaction, have been recently 
studied both experimentally and theoretically. In fact, a crystal-field in- 
teraction effects on the transition temperature are investigated by several 
methods such as effective field theory [H O [Sj Hj , finite cluster approximation 
[S] , mean field theory [6] , Migdal-Kadanof renormalisation group method [7] 
and cluster variational method [8|. However, there are some disagreements 
among those theoretical studies such as the existence of tricritical points 



*benyous@fsr.ac.ma 
^e-mail: elkenz@fsr.ac.ma 



and other features. Experimentally, the MnNi{EDTA) — QH2O complex 
has been shown to be an example of a mixed spin system [9] . Another inter- 
est towards the mixed spin Ising models can be related to the modelling of 
magnetic structures suitable for describing a ferrimagnetism of certain class 
of insulating materials. Indeed, these systems provide simple but interest- 
ing models to study molecular magnetic materials that are considered to 
be possibly useful materials for magneto-optical recordings. Furthermore, 
important advances have been made in the synthesis of two and three di- 
mensional ferrimagnets, such as 2d organo-metallic ferrimagnets [lUt 
2d networks of the mixed-metal materials [P{Phenyl)/i\[MnCr{oxalate)^n 
|12j . However, the possibility of many compensation points in a variety of 
ferrimagnetic systems has been clarified theoretically \i?>\ [Tl] . 
Since the ferrimagnetic order plays an important role in these materials, the 
investigation of ferrimagnetism in mixed spin systems has rapidly become 
a very rich field of research. As it is well known, these systems have less 
translational symmetry than their single-spin counterparts since they consist 
of two interpenetrating and non equivalent sublattices. From a theoretical 
point of view, many different methods have been employed in these studies. 
In particular, the mixed spin— 1/2 and spin— 1 Ising model has been solved 
exactly in special cases \lb\ I16|. On the other hand, the use of approximate 
methods such as mean field theory, free-fermion approximation, effective 
field theory, high-temperature series expansions, renormalisation group and 
Monte Carlo simulation, have revealed interesting results. 
More recently, new magnetic properties and compensation behaviour are 
found for mixed spins in the presence of a crystal- field [T7]. Such magnetic 
compensation behaviours were predicted by Nel theory of ferrimagnetism 
[Is] , but not revealed in previous works. 

On the other hand, the effective-field theory study, and multicritical be- 
haviour, of the mixed spin Ising model with different anisotropics have re- 
vealed contradictions with earlier works obtained within the same theoretical 
framework. But the mean-field theory study based on the Bogoliubov in- 
equality for the Gibbs free energy [19| , it has been shown that in the presence 
of a single-ion anisotropy, the phase diagrams obtained exhibit a variety of 
multicritical points such as tricritical points and isolated critical points. 
On the other hand, some recent experimental studies performed on various 
single-crystal samples RVO^ [R = La, Nd, Sm, Gd, Er, and Y) [20], showed 
a temperature-induced and magnetization reversal. 

The purpose of this work is to study, via mean field approximation (MFA), 
the infiuence of crystal-field disorder on the phase diagrams and magneti- 
zations of a mixed-spin ferrimagnetic Ising system, in which the two inter- 
penetrating sublattices have spins a = ±1/2 and S = ±3/2, ±1/2. The 
most interesting result emerging from this study is the appearance of new 
types of phase diagrams, in the particular case of two- valued distribution 
of crystal-field. Consequently, three topologically different types of phase 
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diagrams occur. This paper is organized as follows: in the section 2, we in- 
troduce the model and give the details of the MFA. The ground-state phase 
diagram is discussed in section 3. In section 4 we present and discuss our 
results. Finally, section 5 is devoted to summarizes and conclusions. 

2 Model and method. 

Since the MFA method neglects correlations between different spins, it is 
interesting to study the behaviour of complex spin systems, such as the 
ferrimagnetic mixed Ising models. The model we are studying consists of 
two interpenetrating sublattices. One sublattice has spins a assumed to take 
the values ±1/2, the other sublattice has spins S that can take four values: 
±3/2, and ±1/2. The spins S have only the spins a as nearest neighbours 
and vice versa. The interaction between the spins a and S is assumed to be 
an antiferromagnetic exchange. The Hamiltonian of this model is written 
as: 

N/2 

n = JY. ^iS, + J2 ^iSf (1) 

<ij> i=l 

where N is the total number of lattice sites. The exchange interaction 
parameter J is assumed to be positive. The first summation is carried out 
only over nearest pairs of spins and Aj is a quenched random crystal field 
distributed according to the probability distribution [211 [221 [23 [211 • 

r{Ai) = ^[6{A, - A(l + a)) + 6{A, - A(l - a))] (2) 
where a is a positive constant. 

An analogous probability distribution has been used to investigate the crit- 
ical behaviour of '^He He mixtures in random media (aerogel) modelled 
by the spin— 1 Blume-Capel model. In this model, the negative crystal- field 
value corresponds to the field at the pore-grain interface and the positive 
one is a bulk field that controls the concentration of ^He atoms [251 1261 [27] . 

The variational principle based on the Gibbs-Bogoliubov inequality for 
the free energy per site is described by [281 l29l [30]: 

<^ = -Tin{Zo)+ <n-no>o (3) 

Let us denote by and hs the molecular fields associated with the order 
parameters rua =< cr >o and ms =< S >o, respectively, expressed as: 

z 

K = JJ2 < >o= zJms (4) 
i=i 



3 



hs = J ^ < cr >o= zJuia- (5) 

i=i 

where z is the number of nearest neighbours and < ... >o= ^^expi-^pno)°^ 
denotes the average value performed over the Hamiltonian Ho- 
The effective Hamiltonian of the system is given by: 

N/2 N/2 N/2 

1=1 i=l i=l 

The partition function generated by the above Hamiltonian is : 
Z„ = Tr{ea;p(-^)) 

where T is absolute temperature and the Boltzmann's constant has been set 
to unity. The total free energy is given by: 

NJzmams NhaTricr Nhsms 
2 2 2 

-NT J Log{Zo)r{Ai)dAi. (8) 

After the integration over the probability distribution, the free energy per 
spin is given by: 

$ NJzrriams Nhania Nhsms t K 
N = 2 2 2 jl^'"*^'"'''^' 

+i(L09(2e.p(-5^?(l±i!i).„,„(^) + 2e.p(-:Sl±5i)c„.M|)) 

i„,(2e.K-5^?fi^)c<,./.(^) + 2exK-^^i^)c<,.M^)))l(9) 

In order to investigate the magnetizations of the system, the order pa- 
rameters and ms are defined by minimizing the free energy. Then, the 
mean-field equations of state are expressed as follows: 

1 z 

rria = --jtanh{—Jms) (10) 
ms = -^[^ + ^] (11) 

where, 

A = 3exp{— — {l + a))smn{ — - — ) + exp[— — {l + a)))smn{- 



4 



B = 2exp{--{l + a))cosh{^) + 2exp{--{l + a)))cosh{^) 

C = 3exp(-^(l-a))sm/i(^^^) + exp(-^(l-a)))sm/i(^^) 

D = 2exp{- — {l-a))cosh{^^) + 2exp{- — {l-a)))cosh{^^) 

In the above equations, and in all the following, t and d denote the reduced 
temperature T/J and the reduced crystal field A/ J, respectively. The co- 
ordination number z is set to be 4 (square lattice). 

The solutions of the Eqs. [TUHTT] are not unique, the stable ones are those 
minimizing the free energy Eq. El while the others are the unstable ones. If 
the order parameters are continuous (discontinuous), the transitions are of 
second (first) order. 



3 Ground state 

The ground state phase diagram of the system under investigation, according 
to the distribution of the crystal field law (Eq. [2]), is illustrated in Fig. 1. 
Indeed, for very low temperatures and depending on the values of a > 
and the reduced crystal field d = A /J. Eqs. [TOlIllI lead to three solutions: 
(m^r = -1/2, ms = 3/2), (m^ = ^,ms = ^), and {m„ = ^,ms = 1). 
By comparing the energies for all possible configurations, we established the 
ground state phase diagram. This phase diagram is drawn in the reduced 
{d, a) plane. One can distinguish four cases: 

1- For a = 0, a first order transition between the phase (rrio- = — 1/2, ms = 
3/2) and the phase {rria = ^,ms = \) occurs at d = +1. 

2- For < a < 1, two first-order transition lines occur: in one hand between 
the phase {rria = — 1/2, ms = 3/2) and the phase {m^ = ^,ms = 1) 
according to the equation line d = z/(4(l -|- a)) and in the other hand 
between the phase {m^ = ^,rns = l)and the phase (mg- = =^,ms = \) 
according to the equation line d = z/(4{l — a)). 

3- For a = 1, a first order transition between the phase (mo- = ^,ms = |) 
and the ferrimagnetic phase (m^ = ^,ms = 1) appears at d = h. 
4- For a > 1, two first-order transition lines between the (mg- = ^,ms = 1) 
phase and the [m^ = ^,ms = |) phase and the (mg- = =^,ms = |) 
phase and the {ma = ^,ms = 1) phase occur at d = 2^/(4(1 — a)) and 
d = z/ (4(1 -|- q)), respectively. 

It is worth to note that for d = 0, only the phase {m^ = ^,ms = |) 
is stable at T = OK. For a higher temperature, the phase diagrams for 
different values of a and d will be discussed in all the following. 
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Fig. 1 
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Figure 1: The ground state phase diagram estabhshed in the (d = A/J, a) 
plane. (—1/2, 1/2), (—1/2, 3/2) and (—1/2, 1) are the only stable phases for 
very low temperatures. 

4 Phase diagrams and discussions 
4.1 Phase diagrams 

A detailed discussion dealing with finite temperature phase diagrams is dis- 
cussed in this section. For this purpose, we solve numerically the Eqs. [TOl [TT] 
and El A rich variety of phase transitions is observed both when varying 
a in the (tc = Tc/J,d = A/ J) plane, and d in the {tc,a) plane. Indeed, 
the critical temperature is plotted as a function of d for a = (Fig. 2a), 
a = 0.5 (Fig. 2b), a = 1 ( Fig. 2c) and a = 2 (Fig. 2d). In Fig. 
2a ,a = 0, the paramagnetic (0, 0) and ferrimagnetic phases are separated 
by a second-order transition line (solid line). For very low temperatures, 
we found a first-order transition line (dashed line) separating the phases 
{ma- = -Y'j'ms = |) and {rria = ^,ms = ^)- This first-order line is termi- 
nated by an isolated critical point located at {dtr = 1, Ur = 0.099). Whereas, 
for a = 0.5 see Fig. 2b, we found two first-order transition lines. In one 
hand, a first-order transition line separating the phases (rrio- = ^,m5 = |) 
and (rrio- = ^,ms = 1), and terminated by an isolated critical point located 
at {dtr = 2/3, ttr = 0.15). On the other hand, a first-order transition line sep- 
arating the phases {rria = ^,ms = 1) and {rria = ^,1113 = |), terminated 
by an isolated critical point located at {dtr = 2,itr- = 0.045). Above these 
points, a continuous passage to the paramagnetic phase occurs. In Fig. 2c, 
plotted for a = 1, the phase (rrio- = ^,ms = \) disappears at low tempera- 
ture and the only first-order transition line, present in this region, separates 
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the phases (m^ = ^,ms = |) and {ma- = ^,ms = 1), and terminated by 
an isolated critical point located at {dtr = 1/2, tfr = 0.034). Therefore, the 
second order transition-line temperatures arc greater than those found for 
a = and a = 0.5. This phenomenon is still present for a higher value of a 
(see Fig. 2d for a = 2). Indeed, the ferrimagnetic phase (m^ = = 1) 

emerges for large and positive values of the crystal-field. This is due to a 
competition between positive and negative values of the crystal-field which 
favours the ferrimagnetic phase. Consequently, the system exhibits two 
first-order transition lines, for very low temperatures, terminated by two 
isolated critical points: {d = —l,t = 0.11) and (d = 1/3, t = 0.19) separat- 
ing the phases {nia = ^,ms = 1) and (rrio- = ^,ms = |), and the phases 
{rria = ^,ms = |) and (m^ = ^,ms = 1) respectively. 
To give more details concerning the first-order transition lines, we have 
developed, at low temperature, the free energy and entropy of the system. 
Indeed, the free energy and entropy are plotted as function of d for t = 0.035 
and two a values: 0.25 and 1, in figures 3a and 3b, respectively. In accor- 
dance with Figs. 26 and 2c, we observe a discontinuous change of the free 
energy slope at first-order transition temperatures. Consequently the en- 
tropy is discontinuous at these temperatures (see inset of Figs. 3a and 3b). 
For the second-order transition lines, we have plotted the total free energy 
and entropy, versus temperature, for a = 0.25 and d = —3 (see Fig. 3c). 
We can remark that the entropy slope varies discontinuously at a second- 
order transition temperature. As consequence, the specific heat will exhibit 
a discontinuity at this transition temperature. 

In order to outline the effect of system behaviour as a function of the 
parameter a we plot, in Figs. 4, the transition temperatures tc for several 
values of the crystal field d. Indeed, for d = 2, Fig. 4a shows the existence 
of a second-order transition line between the paramagnetic and the partly 
ferrimagnetic phases and a first-order transition between the phases {nia- = 
^,ms = 5) and (m^ = ^,ms = 1) terminated by an isolated critical 
point (a = 3/2, t = 0.055). Whereas for d = , see Fig. 4b, the second 
order transition line from the ferrimagnetic phase to the paramagnetic one 
is independent on a values because of the absence of the crystal field. In Fig. 
4c, plotted for d = —1, we found a first order transition line separating the 
phases (mo- = ms = |) and (mg- = ms = 1), which terminates at an 
isolated critical point {a = 2,t = 0.09) between the ferrimagnetic and partly 
ferrimagnetic phases. At high temperatures, a second order transition line 
between these phases and the paramagnetic one is found. 

4.2 Magnetic properties 

In this part, we focus our interest on the magnetizations corresponding to 
the previous phase diagrams as a function of the crystal field d and the 
parameter a, for fixed temperature values. Indeed, Figs. 5a and 5b show 
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the temperature dependence of magnetizations as a function of d for a = 
and d = 0.5 respectively. The former figure presents a first order transition of 
the magnetization ms, for a low temperature t = 0.05: from 3/2 to 1/2 and 
a second order transition of the magnetization 777.5, for a higher temperature 
t = 1.5: from 3/2 to the paramagnetic phase. The second figure, Fig. 5b, 
shows a double first order transition of the magnetization 7775, for a very low 
temperature t = 0.04: from 3/2 to 1 and from 1 to 1/2, and a second order 
transition for a higher temperature t = 1.1: from 3/2 to the paramagnetic 
phase for increasing values of the crystal field d. In Fig. 5c we present 
the behaviour of the magnetization mg as a function of the crystal field for 
a = 1 and two temperature values: t = 0.03 and t = 0.6. In agreement with 
Fig. 2c, we found a discontinuity of this magnetization at a first order point 
between the phases {ma = ^,ms = |) and (777^ = ^,1113 = 1), for a low 
temperature {t = 0.03). On the other hand, a continuous passage occurs 
from the same phases for a higher temperature {t = 0.6). 

The re-entrant behaviour found in Fig. 2d, is well illustrated in Figs. 
6a and 6b plotted, for a = 2, and two temperatures t = 2.2 and t = 2.3, 

respectively. Indeed, Fig. 6a shows that for t = 2.2, the magnetizations 
ma,ms and consequently M = (rria + ms)/2 drop to zero in the region 
< d < 3.5. This is in good agreement with Fig. 2d. On the other hand. 
Fig. 6b shows that for a higher temperature {t = 2.3) this phenomenon is 
inverted in the region —3 < d < 0, so that the magnetizations ma,ms, and 
consequently M = (rria + ms)/2, are nulls out of this region. This is once 
again in agreement with Fig. 2d. 

Let us now show the behaviour of the magnetizations ma,ms, and M = 
{niu + 777s) /2 SiS a function of the parameter a for fixed values of crystal 
field and temperature. For this purpose we plot in Figs. 7a and 7b these 
magnetizations as a function of a for {d = 2,t = 2) and {d = —l,t = 2.5) 
respectively. The first figure shows that for positive values of the crystal 
field, the increasing a values effect is to force the ordered phase, in a region 
of temperatures lower than tc- This is the case in Fig. 7a, for t = 2. For a 
negative value of the crystal field, this phenomenon is inverted as it is shown 
in Fig. 7b, plotted for d = —1. One can note that the results found in Figs. 
7a and 7b, arc in good agreement with those illustrated in Figs. 4a and 4c, 
respectively. On the other hand, it is obvious that the magnetizations are 
independent of the parameter a in absence of the crystal field, as it is well 
illustrated in Fig. 4b. 

To complete this study, we have investigated the effect of increasing the 
exchange interaction parameter J, at very low temperature when keeping a 
constant. Indeed, the results we found for a low temperature t = 0.025 and 
selected values of a, showed that the transitions obtained in the ground state 
phase diagram (Fig. 1) are still present. For a = 0.5, Fig. 8a showed that 
the system can exhibit the phases (777^ = ^,ms = 5), {m-a = -^,1^3 = 1) 
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and {ma = ^,ms = |) when increasing the parameter J at a positive and 
constant crystal field {d > 0). For a = 1.0, the system undergoes a transition 
from the phase (mo- = ^,mg = 1) to the phase (mg- = ^,ms = |) for 
positive and constant crystal field, when increasing the exchange interaction 
J, see Fig. 8b. For a large value of a (a = 2), Fig. 8c shows that the effect 
of increasing the parameter J on the phase transitions (rrio- = ^,ms = 
{rrifj = ^,ms = 1) and (m^ = ^,ms = |) is to displace these transitions 
towards large and positive values of the crystal field. 
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5 Conclusions 



We have investigated a mixed spin £7 = 1/2 and spin 5 = 3/2 Ising model on 
a square lattice, using the mean field approximation {MFA) . The effect of a 
random crystal field on the magnetic properties of the system is investigated. 
Indeed, our results revealed many interesting phenomena, namely, several 
topologically different types of phase diagrams. Furthermore, these phase 
diagrams present rich varieties of phase transitions with first and second 
order phase transition lines. These lines are linked by trier itical points and 
terminated at isolated critical points. Finally, the effect of increasing the 
exchange interaction parameter J, at very low temperature when keeping a 
constant, is investigated. The results we found for a very low temperature 
and sleeted values of a, showed that the transitions obtained in the ground 
state phase diagram are still present, but displaced towards large and posi- 
tive values of the crystal field. 
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Figure 2: The critical temperature as a function of d plotted for: a = 
(a), a = 0.5 (b), a = 1 (c) and a = 2 (d). The full lines correspond to the 
second-order transitions, whereas the dashed lines represent the first-order 
transitions. The tiny circles denote tj® isolated critical points. 



0.0 



-0.2 
-0.4 

-0.6 

S -0.8 

c 

a> 

I -10 
1^ 

-1.2 
-1.4- 

-1.6 



a =0,25 
t=0,035 



-1.6 



a=1,0 





t=0,035 




























0.0 0.5 1.0 1.5 2.0 2.5 



-4.5 




Figure 3: In a and b, the free energy and entropy (from expressions developed 
at low temperature) are plotted versus the reduced crystal-field d for t = 
0.035 and two a values: 0.25 and Ij^gfespectively. In c, the same physical 
quantities are plotted as a function of temperature for a = 0.25 and d = —3. 
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Figure 4: The transition temperatures tc as a function of the parameter a 
for selected values of the crystal field: d = 2 (a) , d = (b), and d = —1 (c). 
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Fig. 4a 
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Figure 5: The dependency of magnetization ms as a function of the reduced 
crystal-file d for a fixed value of a and two reduced temperature values: 
a = and (t = 0.05 and t = 1.5) (a), a = 0.5 and {t = 0.04 and t = 1.1) (b) 
and a = 1. and {t = 0.03 and t = 0.6|5(c), respectively.. 
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Figure 6: The dependency of magnetizations ms, and M = {mfj+rns)/2 
as a function of the reduced crystal- file d for a = 2, and two temperature 
values: t = 2.2 (a) and t = 2.3 (b), respectively. 



16 



0,5 
0,4 
0,3 
0,2- 
0,1 
0,0 
-0,1 
-0,2 



d=2 
t=2 



Fig. 6a 



- m 

- m 

- m -UTi 



,„„0000000000000 



0,5 



2,0 



Fig. 6b 



0,3 



(13 
N 

1) 0,1 
(0 



0,0 





\ 

\ 




d=-1 
t=2.5 


— • — m 

— » — m 
— — m -i-m 

















1,0 1,5 



Figure 7: The behaviour of the magnetizations ma,ms, and M = (m^ + 
771-5 )/2 as a function of the parameter a for fixed values crystal field and 
temperature for {d = 2,t = 2) in (a) and (d = —l,t = 2.5) in (b). 
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Fig. 7a 



12 



10- 


a=0.5 




8- 




(-1/2,1/2) 


6- 

4- 




(-1/2,1) 


2- 
n 




(-1/2,3/2) 



1 2 3 4 5 6 

J 



Fig. 7b 
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Fig. 7c 
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Figure 8: Phase (diagrams in the plane (A, J), at very low temperature t = 
0.025 for Q = 0.5 (a), a = 1 (b) and a = 2 (c). In (a) the system exhibits the 
phases {rria = ^,ms = 5), (m^r = = 1) and {rria = =^,ms = |), in 

(b) the system undergoes a transition from the phase (m^- = ^,ms = 1) to 
the phase {rria = ^,ms = 1), while (c) shows the phases {rrio- = ^,ms = 
{ma = ^,1713 = 1) and {rricr = ^,ms = |) when increasing the 
exchange interaction parameter J for fixed values of the crystal field A. 






